THE EULER-MARUYAMA APPROXIMATION FOR THE 
ABSORPTION TIME OF THE CEV DIFFUSION 
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Abstract. A standard convergence analysis of the simulation schemes for the 
hitting times of diffusions typically requires non-degeneracy of their coefficients 
on the boundary, which excludes the possibility of absorption. In this paper we 
consider the CEV diffusion from the mathematical finance and show how a weakly 
consistent approximation for the absorption time can be constructed, using the 
Euler-Maruyama scheme. 



1. Introduction 

In scientific computations, expectations with respect to probabilities, induced by 
continuous time processes, are often replaced by Monte Carlo averages over inde- 
pendent trajectories. For diff'usions, generated by stochastic differential equations 
(SDEs), the trajectories are usually approximated numerically (see e.g. [E]). The 
accuracy assessment of such numerical procedures is a well studied topic and the 
available theory establishes and quantifies the convergence of the approximations to 
the actual solution in a variety of modes, depending on the properties of the SDE 
coefficients. This in turn typically suffices to claim convergence of expectations for 
path functionals, continuous in an appropriate topology, but unfortunately, may not 
apply to discontinuous functionals, some of which arise naturally in applications. 

One important example of such a functional is the hitting time of a domain 
boundary. Let X = {Xt)t>o be the diffusion process on M, generated by the ltd SDE 

dXt = h{Xt)dt + a{Xt)dBu = x E M, (1.1) 

where B = {Bt)t>o is the Brownian motion and the coefficients b{-) and a(-) are 
functions, assumed to satisfy the regularity conditions, guaranteeing existence of 
the unique strong solution (see e.g. (23j[23]). The hitting time of the level ^ E M is 

Ti{X) := inf{t >0:Xt=e}, 

where inf{0} = oo. Thus Te{X) is an extended random variable with values in the 
Polish space M4. := M+ U {00}, endowed with the metric p{s,t) = | arctan(s) — 
arctan(i)|, s, t E M+. 

Consider a family of continuous processes X^ = {X^)t>o, generated by a numerical 
scheme with the time step parameter 6 > (such as e.g. the Euler-Maruyama 
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recursion ()2.2p and ()2.3p below) and suppose that X approximates the diffusion X 
in the sense of weak convergence. More precisely, let C([0,oo),M) be the space of 
real valued continuous functions on [0, oo), endowed with the metric 

q(u,v) =^ 2~U\iY>\ut - vt\, n,u G C([0,oo),M). (1.2) 

Then X^ converges weakly to X, if for any bounded and continuous functional 
h : C([0,oo),M) ^ M 

lim E/i(X^) = E/i(X). (1.3) 

(5->-0 

Such convergence can often be established using the techniques, developed in e.g. 

0, [n], [m. 

In particular, if (f) : C([0,oo),M) i— )• M_|_ is a functional, almost surely continuous 
with respect to the measure induced by X, then (|1.3p implies 

runEf{^{X'))=Ef{cP{X)) (1.4) 

for any continuous and bounded function / : M_|_ i— )• M. In other words, the weak 
convergence of the processes X^ — >• X implies the weak convergence of the random 
variables (t>{X^) — )• <f>{X) or, equivalently, the convergence of the probability distri- 
bution functions lim5_).o P(<^(X'') < x) = P(0(X) < x) for any x G ]R_|_, at which the 
distribution function P{(j){X) < x) is continuous. 
Let us now take a closer look at the hitting time 

Ti{u) = mf{t >0:ut = i}, n G C([0,oo),M), 

viewed as a functional on C([0, oo), M). Clearly ti{-) is discontinuous at some paths: 
for example, if Ut = max(l —t, 0) and ul = max(l — t, 1/j), then Q{ut,ul) < 1/j — )■ 
as j — )• oo, but r(n) = 1 and t(u^) = oo for all j. 

On the other hand, T£{-) is continuous at any u, which either upcrosses i: 

T(:^{u) := inf{i > : Ut > ^} = T£{u), if uq < i 

or downer osses i: 

T£-{u) := mi{t > : Ut < i} = T£{u), if uq > £. 

If this type of paths is typical for the diffusion (jl.ip . i.e. if the set of all such paths 
is of full measure, induced by the process X, then Ti{-) is essentially continuous and 
the weak convergence X^ — t- X still implies the weak convergence T£{X^) — t- Ti{X). 
However, if i is an accessible absorbing boundary of X, the paths which hit i cannot 
leave it at any further time. For such diffusions Te{-) is discontinuous on a set of 
positive probability and the weak convergence Ti{X^) — )■ T£{X) cannot be directly 
deduced from X^ — > X. 

Hitting times play an important role in applied sciences, such as physics or finance, 
and since their exact probability distribution cannot be found in a closed form 
beyond special cases, practical approximations are of considerable interest. There 
are two principle approaches to compute such approximations. One is based on the 
fact that the expectation of a given function of the hitting time solves the Dirichlet 
boundary problem for an appropriate PDE, and thus the approximations can be 
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computed using the generic tools from the PDE numerics. Sometimes, the particular 
structure of the emerging PDE can be exploited to calculate expectations of special 
functions of hitting times, such as moments (as e.g. in the linear programming 
approach of [TO]). 

The probabilistic alternative is to use the Monte Carlo simulations, in which the 
diffusion paths are approximated by numerical solvers. Typically the diffusions are 
simulated on a discrete grid of points and the evaluation of the hitting times requires 
construction of the continuous paths through an interpolation. The naive approach 
is to use the general purpose interpolation techniques, such as the one used in our 
paper (see ()2.3p below). Better results are obtained if the possibility of having a hit 
between the grid points is taken into account as in e.g. [T7], [6], [7], [T2] . 

The convergence analysis of the approximations of the hitting times, based on the 
various numerical schemes, appeared in [19], [211 [20], [H]) [H]- The results obtained 
in these works, assume ellipticity or hypoellipticity of the diffusion processes under 
consideration, which corresponds to the case of non-absorbing boundary in the pre- 
ceding discussion. The analysis beyond these non-degeneracy conditions appears to 
be a much harder problem. 

In this paper we consider a particular diffusion on M4., with an absorbing bound- 
ary at {0}. As explained above, the absorption time in this case is a genuinely 
discontinuous functional of the diffusion paths, which makes the convergence anal- 
ysis of the approximations a delicate matter. We propose a simple approximation 
procedure, based on the Euler-Maruyama scheme, and prove its weak consistency. 

In the next section we formulate the precise setting of the problem and state the 
main result, whose prove is given in Section [3l The results of numerical simulations 
are gathered in Section [H and some supplementary calculations are moved to the 
appendices. 



2. The setting and the main result 
Consider the diffusion process X = {Xt)t>Q, generated by the Ito SDE 

dXt = i-iXtdt + aXfdBt, x e M+, t>0, (2.1) 

where /i € M, cr > and p G [1/2, 1) are constants and B = {Bt)t>o is the Brownian 
motion, defined on a filtered probability space (ri, 3", (3"j), P), satisfying the usual 
conditions. This SDE has the unique strong solution (Proposition 1 in [1]) and is 
known in mathematical finance as the Constant Elasticity of Variance (CEV) model 
(see e.g. [3]). For p = 1/2, it is also Feller's branching diffusion, being the weak 
limit of the Galton - Watson branching processes under appropriate scaling. 

The process X is a regular diffusion on R_|_ U {0} and a standard calculation 
reveals that {0} is an absorbing (or Feller's exit) boundary (see §6, Ch. 15 in |13j). 
We will denote by (Px)a;GK+ the corresponding family of Markov probabilities with 
PxiXo = x) = 1, induced by X on the measurable space C([0, oo),M+) with the 
metric (|1.2p . 
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Consider the continuous time process X = {Xf), t > 0, which satisfies the Euler- 
Maruyama recursion at the grid points tj G (5^4- 

4=4-i+^K-i)^'^ + ^K-i)^''^^-^' Xo = x>0 (2.2) 
and is piecewise hnear otherwise: 

^f = ^^4 + ^4-i' *e (2-3) 

where x"*" := max(0,x), 5 > is a small time step parameter and (^j)j6z+ is a 
sequence of i.i.d. iV(0, 1) random variables. 

Since the diffusion coefficient of ()2.ip degenerates and is not Lipschitz on the 
boundary {0}, this SDE does not quite fit the standard numerical frameworks such 
as [15] or [22]. Nevertheless the scheme (j2.2p does approximate the solution of 
()2.ip in the sense of the weak convergence of measures, as was recently shown in 
[1] (see also [27], [26], [2]). Consequently, for any Pa;-a.s. continuous functional 
(p: C([0,oo),M) h^M+U{oo} 

0(X^)f->0(X), (2.4) 

where ^ stands for the the weak convergence, defined in ()1.4p . Since a typical 
trajectory of X oscillates around the level a > after crossing it, the functional Ta{-) 
is Pi;-a.s. continuous for a > (see Lemma 13.31 below) and hence Ta{X^) ^ Ta{X) 
as — > 0. 

This argument, however, does not apply to ro(-), since it is essentially discon- 
tinuous, as discussed in the Introduction. Leaving the question of convergence 
To{X^) ^ To{X) open, we shall prove the following result 

Theorem 2.1. For any /3 G ^0, Y~f)? 

Ts,{X')^ToiX), 6^0. (2.5) 

Note that T^fi{-) is a continuous functional for any fixed 5 > and hence can be 
seen as a mollified version of the discontinuous tq{-). The parameter /3 controls the 
mollification, relatively to the step-size parameter of the Euler-Maruyama algorithm. 

In practical terms, the convergence (j2.5p provides theoretical justification for the 
procedure, in which the approximate trajectory X^ , generated by (j2.2p and (|2.3p . is 
stopped not at 0, but at 6^ > 0, which only approaches zero as (5 — )■ 0. Our method 
does not quantify the convergence in terms of e.g. rates, but the numerical exper- 
iments in Section [5] indicate that this stopping rule produces practically adequate 
results. 

Remark 2.2. As will become clear from the proof, our approach exploits the local 
behavior of the SDE coefficients near the boundary and can be applied to the more 
general one-dimensional diffusions of the form (jl.ip , for which a weakly convergent 
numerical scheme is available. For example, all the arguments in the proof of Theo- 
rem [2T] directly apply to the diffusion, whose coefficients b{x) and a{x) have a similar 
local asymptotic at as the CEV model. The SDE (12. ip is a study case, which seems 



5 



to capture the essential difficulties of the problem, related to the degeneracy of the 
SDE coefficients on the absorbing boundary. It is a convenient choice, since the 
weak convergence ()2.4p . being the starting point of our approach, has been already 
established for the CEV model in [T], and the explicit formulas for the probability 
density of to{X) are available, allowing to carry out the numerical experiments. 

3. The proof of Theorem 12.11 

The proof, inspired by the approach of S.Ethier is based on the following 
observation (a variation of Billingsley's lemma, see Proposition 6.1 |3]) 

Proposition 3.1. Let Q", n G N and Q be Borel probability measures on a metric 
space S such that Q" A Q as n —7- oo. Let S' be a separable metric space with 
metric p, and suppose that (j))., k £ N and (p are Borel measurable mappings of S 
into S' such that 

(i) is Q-a.s. continuous for all k £ N 

(ii) limfc (pk = 4>, Q-a.s. 

(iii) for every rj > 0, 

limIh^^Q"fp(0fe,(/.„J > r/) = 0. 

k n \ J 

for an increasing real sequence (a„). 
Then Q"^ o (f)~^ ^ Q o (f)~^ as n —t- oo. 

Proof. Let hhe a continuous with respect to p bounded real valued function on S', 
then 

\E^h{cP^J - Eh{(P)\ < E^\hi(baJ - h{(bk)\ + 

\E^h{^k)-Eh{(l)k)\ + E\h{M-h{(P)\, 

where E and E" denote the expectations with respect to Q and Q" respectively. 
Since h is continuous and (pk — t- (p, Q-a.s., the last term vanishes as /c — t- oo by the 
dominated convergence. Moreover, since for any fixed k, (pk is Q-a.s. continuous and 
Q" ^ Q as n — )• oo, 

lim|E"/i((/.fc) - Eh{(t)k)\ = 0, V/c. 

Consequently, 

IE^IE^|E"/i((/)«J - Eh{(p)\ < ]hR]h^E'^\h{(j)a„) - h{(pk) \ = 

k n ' ' k n ' ' 

where the latter equality holds by (jni|) . The claim follows by arbitrariness of h. □ 

Let us now outline the plan of the proof. In our context, the probability measures, 
induced by the family , play the role of Q'' and by Proposition 13.21 they converge 
weakly to the law Q := of the diffusion X. Since the diffusion coefficient of 
()2.ip is positive, away from the boundary point 0, Tir{-) s. continuous 

functional (Lemma 13. 3p and hence ^ of Proposition [STTl holds. On the other hand, 
lime_^o "^e (^) = ''"o(^^) for any u G C([0,oo),M) (Lemma 13. 4j) . which implies ([n]) of 
Proposition 13.11 The more intricate part is the convergence (Iml) , which is verified 
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in Lemma 13.51 using the particular structure of the SDE p. II) . The statement of 
Theorem 12.11 then follows from Proposition 13. 1[ 

The following result is essentially proved in [l]: 

Proposition 3.2. The processes (X^), defined by ([22]) and ([23]), converge weakly 
to the diffusion X, defined by (12. ip . as 5 — > 0. 

Proof. For the process, obtained by piecewise constant interpolation of the points 
generated by the recursion (12. 2p . the claim is established in Theorem 1.1 in [T]. The 
extension to the piecewise linear interpolation (12. 3p is straightforward. □ 

Lemma 3.3 (Proposition 4.2. in j4]). For allx > e > 0, ''"e(-) is a Px-a.s. continuous 
functional on C([0,oo),M). 

Proof. We shall prove the claim for completeness and the reader's convenience. For 
u G C([0, oo), M) and £ > 0, let 

T£-{u) := inf{t > : lif < e} 

and define 

A:={u£ C([0,oo),M) : T^(n) = t^-(u)} . 
We shall first show that u i— t- r^iu) is continuous on A, i.e. that 

£>K,n)^0 =^ /)(r,(n"),Te(u)) ^0, e A, (3.1) 
and then check that 

Px{XeA) = l \fx>e. (3.2) 

To this end, note that if u G ^ and Te(ti) = oo, then m.mt<T ut > £ for all T > 
(recall uq > e for u £ A). Thus g{u'^,u) — )• implies lim„ mint<7^ > e and hence 
lim „ Tpju^) > T. Since T is arbitrary, lim„r£(ti") = oo, i.e. ()3.ip holds. 

Now take u £ A, such that T^iu) < oo. If T£(u) = 0, the claim obviously holds by 
continuity of u. If r^iu) > 0, fix an 77 > such that r^iu) — r] > 0. Since 

mm Ut > £ and sup \Uf — Ut\ > U 

t<Te{u)-ri t<Te(u)-ri 

we have lim„ minj<^^(„)_^ > e and thus lim „ T^ju^) > r^iu) — rj. On the other 
hand, as Te{u) = Te_(u), for any i] > 0, there is an r < Ts{u) + t], such that Ur < e. 
It follows that lim„ M„,(r) < e and hence lim.„ T^iuy^ < Te{u) + rj. By arbitrariness of 
r], we conclude that lim„ Te(u„) = Te(M) and (j3.ip holds. 

It is left to show that (j3.2p holds. The diffusion X satisfies the strong Markov 
property and thus (we write for Ts{X)), 

Px{A) = P^r, = T,_) = E,1{,^<^}P,(t,_ = 0) + P,(r, = 00), 

and hence the required claim holds, if P£(t£„ = 0) = 1. 

Take now / to be a scale function of the diffusion X, i.e. a solution to the equation 
Lf = 0, where L is the generator of X: 

{Li)) (x) = nxip'ix) + ^a'^x'^Pij"{x), X > (3.3) 
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It is well known, e.g. [H], that we can take it to be positive and increasing, specifi- 
cally, for p.3p . 



can be taken. The process f{Xt) is a nonnegative local martingale and thus a 
supermartingale (e.g. [H] p. 197). The random variable t A r^- is a bounded 
stopping time and by the optional stopping theorem we have for any t > 

Ee/(XiA.,_)</(e). (3.4) 

By the definition of Te~ and path continuity of X, it follows that Xt^r,^ ^ £) and 
since / is increasing, we have that f{Xt/\T^_) > f{e). Thus it follows from (|3.4|) that 
Pf;{Xt/\r^_ = s) = 1 for all t > and, consequently, [X, X]t/\r^_ = 0, P^-a.s. for 
t > 0. On the other hand, [X,X]r^_ = a^Xs^ds > 0, on the set {ts- > O}, 
Pg-a.s. The obtained contradiction implies Pe(Te_ > 0) = 0, as claimed. 

□ 

Lemma 3.4. lim^^Q Ti:{u) = tq{u) for any u € C([0, oo), R+). 

Proof. If To(ti) = oo, then Ta.m.t<T ut > ^ and hence lim^ ^^g Tg^^u) > T and the claim 
follows, since T is arbitrary. If tq{u) < oo, then for an ?] > 0, inf4<^y(„)_^ > and 
thus lim^ ^Q Tp:(u) > tq{u) — Tj. For sufficiently small e > 0, Tfr{u) < tq{u) and the 
claim follows. □ 

Let P be the probability on the space, carrying the sequence (^j) (see (|2.2p ) and 
denote by (Px)a;>o the Markov family of probabilities, corresponding to the discrete 
time process {Xf.) with ¥j.{Xq = x) = 1. Since the process (X^) is piecewise linear 
off the grid the condition (jui]) of Proposition 13.11 follows from 

Lemma 3.5. For any j3 G ^0, and i] > 

limm¥jp{T,iX'),Ts,iX')) > 7?) = 0. (3.5) 

Proof. Roughly speaking, (j3.5p means that a trajectory of X^, which approaches the 
boundary {0}, is very likely to hit it. This seemingly plausible statement is not at all 
obvious, since the coefficients of our diffusion decrease to zero near this boundary, 
making it hard to reach. By letting the level 6^ decrease to zero at a particular rate 
allows to approximate expectations of the hitting times t^i3{X^) by those of Tjti{X), 
which in turn can be estimated using their relations to the corresponding boundary 
value problems. 

In what follows, C, Ci,C2 etc. denote unspecified constants, independent of e and 
5, which may be different in each appearance. Define the crossing times of level a 



^a = 5 inf{j>l:aG (X(_^,X(,]}, aG 



with inf{0} = oo. The sequence (Xl ), j G is a strong Markov process and Ua/^ 



is a stopping time with respect to its natural filtration. 
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Since X is piecewise linear, 

\t,{X^) - Ue\ < 6 and \Tsfi{X^) - i^s^ \ < 6 on the set \^Te{X^) < 00} 
and thus by the triangle inequality 

p{t,{X'),ts,{X')) < p{u„us,) +26. 
Since Te{X^) < Tg^iX^) for 2; > e > 5^, it follows 

¥,(^p{Te{X^),Ts,{X^)) > v) =r,[p{Te{X^),Ts0{X^)) > vM^^) < Oo) < 

T^s^) > 1]- 26, z^e < 00) = E^l{^^<^}P^i (p{0, v^fi) > 1]- 25] < 



lim lim sup E^l u^fj A z^i ) =0, (3-6) 



sup {i^sP > r]') , 

2G[0,£] 

where r]' := tan(?7/2) (assuming 6 is small enough). Further, 

F,{us.B > ri') = ¥,[usp > ?/, z/i > r/') + P^z/^^ > ?/, z^i < 7?') < 

and thus (13.50 holds, if we show 



and 

lim lim sup Pzfi^^s > ^^i ) =0. (3.7) 

Proof of (|3.6p . We shall use the regularity properties of the function 

^(x) := E,To(X) Ari(X), 2; G [0, 1] 

near the boundary point 0, summarized in the Appendix |X1 In particular, ip is 
continuous on the interval [0,1] and is smooth on (0,1]. Note, however, that the 
derivatives of tp explode at the boundary point 0, which is related to the possibility 
of absorption. We shall extend the domain of ip to the whole M by continuity, setting 
ip{x) = for x G ]R\ [0,1]. 

Consider the Taylor expansion 

n 



n 



2 



n 



3! 
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where is between and X^.. After rearranging terms, the latter reads 

n 

where L is the generator defined in l\'A.'6\i . the second term is the martingale 




and the last term is the residual 

Consequently, for an integer k>l and the stopping time v := t'^/s A z/i, 

(i//M/t)-l 

r(5), (3.8) 

where 

r(<5) := V(^^Afe5) - V'(^(tAfc5)-5) - (^./Mfe - M(,/Mfe)_i), 

is the residual term, which accommodates the possible overshoot at the terminal 
crossing time v. 

Recall that Lip = —1 for x G (0, 1) and hence 

i//5Afc-l 

Y (£V)K.J5=-(^^AM)+5. (3.9) 

By Lemma |A.2|, sup^g^Q (^xP\ip'{x)\ + x'^^l'ip" {x)\j < C and hence the increments 
of M„ satisfy 

E^(^|Mj -Mj_i||j|_i) < C < oo, on {j < z>/(5A fc} 

and by the optional stopping theorem (Theorem 2 of §2, Ch. VII, [25]) E^Mj^/^^ji. = 
for z e [0,1]. 

Now we shall bound the residual terms in p.Sp . By Lemma I A. 21 

sup x'^\iIj" {x)\ < oo 
xe(o,i) 
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and hence 



(1) 



u/5Ak 



<E,Y1 l\r{X^^_,)W{x!^j's^ <CidE,iu Ak5). (3.10) 



Similarly, sup^.g(o,i) x^~^^\Tp" {x)\ < oo and by Corollary IB. 21 (applied with r = 0) 



(2) 

i//Mfc-l 



u/5Ak 

ms^/^ <C2V6(E^{uAk6) + 6y (3.11) 



To bound l-R^y^/^^.i | > note that by Lemma E21 



r{x')\<cs\xu 



-2p-l 



and thus 



(3) 

u/SAk' 



j<!^/(5Afc 



< 



4 E + k(4..)%-^r) 



< 



j<i^/5Ak 



j<u/SAk 

where the latter inequality holds, since Xf._^ < 1 on the set {j < v/5 /\k}. Since 
Xf is between Xf and Xf 

K_j"'^"'<K-.)"''"'v(4)"^-\ 

For j < we have Xl_^ > and thus on the set {Xl_^ < Xf ], 



(3.13) 



XLJ"^'"^^ <r^('-^). (3.14) 



On the set {Xl_^ > X^] we have 



5 ^3p 



(x(_^(i + M) + <t(x(_J%V^' 



2p+l • 



(3.15) 
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Note that for a; > 0, a > and b, such that xa + x^b > 0, 

x^P _ (xi-P)'P _ 1 {x^-Pg + b-h)^" 



Q,2p (a;i-Pa + a2p \x'^-Pa + b {x^-Pa + b)^P+^ J ' 

Applying this inequahty to (|3.15p on the set < (JK^^^'P)" 2) | we get 



C5 ^ c^W^.s^'l^p ^ 



2^p(/3(l-p)+|) 



C5 C5C - 

+ • - ' 1 \ X 2p+l - 



(^i + 5P{i/2-/3(i-rt)^ < c^5~^^^~P\ 

where the inequahties hold for all sufficiently small 5 > and we used the bounds 
^t,_i > and /3(1 -p)< 1/2. 

Consequently, on the set {X^.__^ > Xf.^ 

J'^ < C-r6~^^'-P^ + ^'^(^p+Di ji^^i > ^|(/3(i-p)-|) I (3_ig) 

Plugging the bounds p.l6p and (|3.14p into ()3.12p and applying the Corollary IB. 21 
we obtain the estimate 



j<u/5Ak ^ J ^ 

2CsE,iu A M)E, (^(^T + {161 > S-h^ ^ (iCjf + 1) + Cg6^/\ 



< 



where 7 = 1/2 — /3(1 — p) > 0. Using the Gaussian tail estimate 
L{|6l>a} 

we get 



E|6l'"^l{|6l>a} < V(4m-l)!!^e-i'^', m > 1, a > 0, 



r(3) 

-^I'/Mfc-l 



which along with (|3.1Up and (|3.1ip yields the bound 

^z\Ru/SAk-i\ < CuE,{u A k6)6"' + Cii6^l\ (3.17) 
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Finally, by Corollary [B21 
and hence by Lemma lA.21 
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Mu/SAk - 



and consequently 

E,r{6)<^l^{xt^ks)+Cl35'^^(l+E,{uAk6)y (3.18) 
Plugging the estimates (j3.9p . (|3.17p and (j3.18p into (|3.8p . we obtain the bound 
E^z^ A k6) (l - C6"i - C5^'^) < Cb^l^ + ^{z) + sup < oo. 

By the monotone convergence, the latter implies 

^ ^ g^^/^ + i)[z) + SUp^g[o^g;S] 

and by continuity of "0, 

lim sup E^i/ < sup i)(z) iztli. 
^^0 2e[0,e] 2G[0,e] 



verifying (|3.6p . 

Proof of (|3.7p . Consider the function (^(x) := Pa:(To > n), whose domain we extend 
to the whole real line by continuity, setting ip{x) = for x < and (p{x) = 1 for 
X > 1. The process ip{Xj), j > satisfies the decomposition (j3.8p . with ijj replaced 
by (p. Taking into account that {Lip){x) = for x G [0, 1] and the bounds from the 
Lemma [A. H a calculation, similar to the one in the preceding subsection, shows that 

E,(^(Xf) - ^{z) < E,((^(xf) - vixt-s)) + C6'/^ + {CS^ + C,51/^)e,z.. 
On the other hand, 

M^sf^ > ^i) (v'(l) - AXuj) + E,^{X,J > P.(z.5, > ui) (l - sup ifiy) 
and thus, for sufficiently small 5 > 0, 



sup P^(z^5/3 > z^i) < C 1 - sup ip{y) X 

z&[0,e] \ y<5f> J 

( sup ip{z)+ sup EJip{xi)-^{xi_s))+5^''^/\l+ sup E, 
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Since </? is continuous on [0, 1] and ip{0) = 0, 



lim sup > lyi) < C inf ip{z) > 0, 



which verifies (|3.7p . 

□ 



Remark 3.6. The condition /3 < originates in the estimate (j3.14p . which is 

plugged into (|3.12|) . The principle difficulty is that the term ^{^tj^i)^^ 
cannot be effectively controlled for greater values of /3 as (5 —t- 0. For example, it is 
not clear how to bound the right hand side of (13.12P already for (3 = 1 and p = 1/2. 



4. Numerical experiments 

In insurance, one is often interested in calculating the probability of ruin by a 
particular time t > 0. For the diffusion (|2.ip . this probability can be found explicitly 
and for p = 1/2, it has a particularly simple form: 

2x 1 



exp 2 7 ' = 

p4Mx)<t) = { V ^^ty (4.1) 

Note that tq{X) has an atom at {+00} when /i > 0: 

P,(ro(X) < 00) = e-2-M/-' < 1. 

Figure[T]depicts the results of the Monte Carlo simulation, in which the probability 
of absorption (j4.ip has been estimated for particular values of the model parameters, 
using M = 10'' i.i.d trajectories, generated by the Euler-Maruyama algorithm (|2.2p 
and (j2.3p . The relative estimation errors 

Err := ^ ^-^-^ \^ ^ ~ ^ x 100% 

PxiTo{X)<t) 

are plotted versus 6 (in the log scale), along with the 99% confidence intervals, 
based on the CLT approximation. The results appear to be practically adequate: 
for example, the accuracy of 0.1% is obtained already with S = 10~^. The positive 
bias of the error is not surprising, since the earlier absorption is more probable for 
the larger threshold 6^ . 

The simulation results also indicate in favor of the convergence 

lim P^ToiX^) <T) = P^ToiX) < T), 



which remains a plausible conjecture. 
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Figure 1. CEV model with ^ = 0, a = I, p = 1/2, x = 1 and 
t = 5. The corresponding exact value of the absorption probability 
is P^(to(X) <t) = 0.6703... 



Appendix A. Some properties of the CEV diffusion 

In this section we summarize some asymptotic estimates of the absorption times 
for the diffusion (j2.ip near the boundary {0}, which are used in the proof of Lemma 
Recall that 



ip{x) := Px{to > n), and ip{x) := E^tq V n 
are the solutions of the following problems respectively (see e.g. |13|): 

(£(^)(x) = 0, xG(0,l), (^(0) = 0, 99(1) = 1 (A.l) 

and 

(£^)(x) = -l, xe(0,l), ^(0)=^(l) = 0, (A.2) 



where L is the operator, given by (13. Sp . The solutions in the class of continuous 
functions on [0, 1], which are twice differentiable on (0, 1) are given by the formulas 

m- 

V^(x) = §11^11^, XG[0,1], (A.3) 
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and 



^{x)=^{x)j [S{1) - S{y))m{y)dy+ 

2{l-ip{x)) r {s{y)-S{Q))m{y)dy, (A.4) 



where 

s{y) = exp ^— j -^z^~'^^dz^ (scale density) 

ry 

S{y) — S{x) = / s{z)dz (scale measure) 

J X 

'^(y) = 2 2!, / \ (speed density). 

a^y^Ps{y) 

Lemma A.l. The function ip is smooth on (0, 1), and 

sup (\v'{x)\ + x^P^^^"{x)\ + x'^P\^"'{x)\) < 00. 
xe(o,i) ^ ^ 

Proof. The scale density s{x) = exp (^~^3(|z2^^^~^^) ^® smooth on (0, 1), and 
sup \(p'{x)\ = ^ sup 5(3;) < 00. 

a;e(0,l) -^^ij - ^(^Uj a,g[o^l] 

From (lA.ip we get ip"(x) = — and thus 

2p-i I I// \ \ ^ \ n \\ ^ 

sup X ^ \'-P [X)\< 2~ I'/' (^J| < 

xG(0,l) a:G(0,l) 

Differentiating ()A.ip . we get 

11^' {x) + /ix(^"(x) + aV^''"V"(^) + ^f^^a;2P(p'"(x) = 0, (A.5) 
and consequently 

sup x^^|(^"'(x)| < 00. 



xe(o,i) 



Lemma A. 2. The function if) is smooth on (0, 1) and 

sup (x2p-i|V'(x)| + x^P\il}" {x)\ +x2P+i|V'"(x)|) < 00, /orpG (1/2,1) 
a;G(0,l) ^ ^ 

and 

( , ^^,M '{x)\+x^'P\i^"{x)\+x^P^^\i;"'{x)\) <oo, /orp=l/2. 

Vlogll/a; VI' ' ' ' ' V 



□ 



sup 

xG(0,l) 
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Proof. We have 

-iP'ix) = ip'ix) I (5(1) - Siy)) m{y)dy - ^{x) - S{x)) m{x) 



-^'{x) [s{y)-S{Q))m{y)dy+{l-^{x))[s{x)-S{G))m{x), 
and, since S{x) is increasing and 5(2;) < 1, 

-\'4!\x)\<\ip'{x)\ j m{y)dy + ^^^^ _ ^^^^ xm{x) 

/•X 

+ \ip'{x)\ / ym{y)dy + xm{x). 
Jo 

For p = 1/2, it follows that 

xG(o,i) ^og{l/x) V 1 



and for p G (1/2, 1) 



Now (|A.2j) implies 



sup x'^P ^\i^'{x)\ < 00. 
xe(o,i) 



and hence for p £ [1/2, 1), 



sup x'^^\^l}"{x)\ < 00. 

xG(0,l) 



Further, differentiating (|A.2p we see that -0 satisfies (jA.SP as well and hence 

sup x^P+^\ij"{x)\ < 00, 

a;G(0,l) 

which verifies the claim. □ 



Appendix B. An inequality 

Lemma B.l. Let (?7j)jeN be a sequence of random variables and N be an integer 
valued random variable. Then for constants a > and c > and an integer k 

k 

^mhk\ < cE(Af A A;)" + ^ E|r?j|l{|^^|>cja}. 

i=o 



Proof. For a fixed integer k > 1 and a real number a > 
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k-l 

ElwAfel = |??j|l{Ar=j} + E|^?fc|l{^f>fc} < 

P(iV = j) + ck'^PiN > k) + Ej2\Vjmr,\>cj'^} + E|r?fc|l||,,|>,fc.} < 

j=0 j=0 
k 

j=0 

n 

Corollary B.2. Let (Cj)jgN an i.i.d. iV(0, 1) sequence and N be an integer valued 
random variable. Then for any a > 0, p > and integer k 

E|Wr < E(7VAA:)" + C„,p, (B.l) 

with a constant Ca,p depending only on a and p. Moreover, for any r > 0, the sum 
= E"=o l^iPl{|?il>r} satisfies 

ESNAk < 2E|eiri||5^|>,}E(iV Ak) + Cp, (B.2) 

with a constant Cp, depending only on p. 

Proof. Applying Lemma iB.ll with c := 1 and r]j := we obtain the inequality 

(|RT]) with 



Further, define Cj '■= ICjl^l{|Ci|>r} ~ E|^i|^l{|^j|>r.}, so that Ed = and Ci's are 
i.i.d. and E|(^ip™ < oo for any m > 1. It follows that 

E^yE^*) <Cm/r, (B.3) 

with a constant Cm- 



18 



P. CHIGANSKY AND F.C. KLEBANER 



Applying Lemma fB. II with a := 1, c := 2E|^i|'P1||^j|>j.} and ijj := Sj, we get 

k 
3=0 



2E|eiri{|5,|>.}E(iVAA;) + ^yi52 p ["l ^ Ci > c/2') < 



2E|ei|l{|6|>r.}E(^^AA:)+Cj]A. E ij^C 
2E|6|l{|ai>r}E(iVAfc) + ^^^r^/l 
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< 



i=o 



□ 
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